PROCEEDINGS 

OF SCIENCE 



QCD thermodynamics with effective models 



OS 

o 
o 

(N 

Or 
<D 

(N 

S 1 

a; 

in: 

(N 
> 

oo 

(N 

q 

as 
o 
as 
o 



X 



Bernd-Jochen Schaefer 

Institut filr Physik, Karl-Franzens-Universitat, A-8010 Graz, Austria 
E-mail: |bernd- jochen . schaef er@uni-graz . at 



Mathias Wagner 

Institut fur Kernphysik, TU Darmstadt, D-64289 Darmstadt, Germany and 

ExtreMe Matter Institute EMMI, GSI Helmholtzzentrumfilr Schwerionenforschung GmbH, 

D-64291 Darmstadt, Germany 

E-mail: jnathias . wagner@physik . tu-darmstadt . de 



Jochen Warn bach 

Institut filr Kernphysik, TU Darmstadt, D-64289 Darmstadt, Germany and 
Gesellschaft filr Schwerionenforschung GSI, D-64291 Darmstadt, Germany 
E-mail: 



jochen . wambach@physik . tu-darmstadt . de 



In this talk we extend the Polyakov-quark-meson model to Nf = 2 + 1 quark flavors and study 
its bulk thermodynamics at finite temperatures in mean-field approximation. Three different 
Polyakov-loop potentials are considered. Our findings are confronted to recent QCD lattice sim- 
ulations of the RBC-Bielefeld and HotQCD collaborations. Furthermore, the finite chemical po- 
tential expansion of the quark-number susceptibility in a Taylor series around vanishing chemical 
potential is analyzed. By means of a novel algorithmic differentiation technique, we have calcu- 
lated Taylor coefficients up to 24 th order in the model for the first time. This allows the systematic 
study of convergence properties of the Taylor series. 
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1. Introduction 

The deeper understanding of strongly interacting matter under extreme conditions plays a 
crucial role in many research programs. In particular, the search for possible (tri)critical endpoints 
in the phase diagram of strongly interacting matter is a major focus of, e.g., the planned CBM 
experiment at the FAIR facility. The underlying theory, QCD predicts several phases which are 
linked to certain phase transitions. For example, the confined, hadronic phase at low temperatures 
and chemical potentials (net quark densities) is separated from the deconfmed quark-gluon plasma 
phase at high temperature and densities. Chiral and deconfinement aspects of these QCD phase 
transitions are the major focus of this talk. 

For small temperatures and finite chemical potentials the chiral phase transition is probably of 
first-order while a crossover is expected at high temperatures and small chemical potentials. This 
suggests the existence of at least one critical endpoint (CEP) where the first-order transition line in 
the phase diagram terminates. At the CEP the transition is of second-order. The location and even 
the possible existence of additional endpoints in the phase diagram is still an open question (see 
e.g. [[!]]). Furthermore, the coincidence of the chiral and the deconfinement transition at vanishing 
chemical potentials and a possible confined but chirally symmetric phase, the quarkyonic phase, 
[^] at finite chemical potentials are much under debate. 

Different regimes of the QCD phase diagram can be explored by employing various theoretical 
methods. Lattice QCD simulations are applicable at zero or imaginary chemical potentials but at 
finite real chemical potentials the fermion sign problem is still a considerable obstacle. At finite 
temperatures recent lattice simulations describe the QCD thermodynamics very trustfully since 
larger volumes and quark masses closer to their physical values are used. Despite this progress 
some lattice groups differ in their predictions even at vanishing chemical potentials. The RBC- 
Bielefeld group sees a coincidence of the chiral and deconfinement transition while the Wuppertal 
group obtains a larger critical temperature for the deconfinement transition [|3|, [l], |5]]. Several ex- 
trapolations techniques towards small finite chemical potentials such as the reweighting method, 
imaginary chemical potential or a Taylor expansion around vanishing chemical potentials have 
been proposed (see ^ for an overview). 

Another theoretical method is based on NJL-type or quark-meson models which incorporate 
chiral symmetry breaking in an effective manner. These models do not suffer from any limitations 
at finite chemical potentials or finite volumes but cannot address the deconfinement transition due 
to the lack of confinement. However, such types of models can be augmented with the Polyakov 
loop resulting in Polyakov-loop NJL-type (PNJL) or Polyakov-loop quark-meson (PQM) models. 
With these extended models both phase transitions become accessible. 

In this talk we report on recent results from a chiral PQM model with Nf = 2+1 quark flavors. 
We compare the bulk thermodynamics of this model for various Polyakov-loop potentials with 
newer QCD lattice data where the larger quark masses on the lattice are also explicitly considered 
in the model comparison. 

Based on a novel numerical method, the algorithmic differentiation technique, higher deriva- 
tives of the thermodynamic potential can be extremely precise calculated [{J. This method, applied 
to the PQM model, allows us to investigate convergence properties of the Taylor expansion method 
used on the lattice. 
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2. Polyakov-quark-meson model for three quark flavors 

The Polyakov-quark-meson model for three quark flavors is a combination of the chiral linear 
a-model with the Polyakov loop <&(x), the thermal expectation value of a color traced Wilson loop 
in temporal direction [Hj. Since the Polyakov loop in these types of effective models is used as a 
classical variable, implementation details of this quantity are not important. It serves as an order 
parameter of the center symmetry in the limit of infinitely heavy quarks. It is finite at high temper- 
atures corresponding to the deconfined plasma phase, where the center symmetry is spontaneously 
broken and it vanishes in the confined, center-symmetric phase. However, in a system with dynami- 
cal quarks the center symmetry is always broken explicitly but <I> still seems to be a useful indicator 
of the confmement/deconfinement transition. The PQM Lagrangian consists of a quark-meson part 
and a Polyakov-loop potential % (<!>,<!>), which depends on the Polyakov-loop variable <I> and its 
hermitian conjugate <!>. The coupling of the background gauge field to the quarks is achieved by 
replacing the standard derivative in the quark-meson contribution with a covariant derivative 
Dfi = dp — iAjj with Afi = 5^oA°. This leads to the Lagrangian 

J? PQM = q(ip-g<t> 5 )q + £' m -W(<5>,<5>) , (2-1) 

where the interaction between the scalar (a a ) and the pseudoscalar (n a ) meson nonets and the three 
quark flavors q is implemented by an Yukawa-type vertex 

A 

= -f (Oa + iYsKa) • (2-2) 

The remaining, purely mesonic contribution reads 

J? m = Tr(d^d^) -m 2 Tr(^^)-X l [Tr(<^)] 2 - A 2 Tr (f r 0) 2 

+c(det(0)+det(0 t ))+Tr[#(0 + t )] , (2.3) 

with the short-hand notation <p = X a (o a + ilt a ) /2. The Polyakov loop potential % preserves the 
center symmetry of the pure Yang-Mills (YM) theory and several explicit choices are possible 
which we compare with each other in the following. 



2.1 Different choices of the Polyakov loop potential 



Several parameterizations of the Polaykov loop potential have been proposed in the last years, 
see Wh. The simplest choice is based on a Ginzburg-Landau ansatz and results in an expansion in 
terms of the order parameter 



poly 



y4 



j 2 (|<J>[ 2 + j<!>| 2 ) -b 3 (& +<l> 3 ) +b 4 (|<£| 2 + |<I>| 2 ) 



(2.4) 



with a temperature-dependent coefficient b2{T) = ao + a\ (^) + a 2 (y) + a 3 ( y) . The parame- 
ters are adjusted to the pure gauge lattice data. An improved ansatz, motivated by the SU (3) Haar 
measure, results in 



log 



^a(T)®<i> + b{T)ln 



1 - 600 + 4 ( 4> J + 0> J ) - 3 ( ) 



(2.5) 
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54e a/T <f>3> + In (1 - 6<M> - 3(<M>) 2 + 4(<I> 3 + <l> 3 )) 



(2.6) 



with the temperature-dependent pre-f actors a(T) =ciq + ci\ (y) ^^(t) an( * b(T) = b$ (jl-) . In 
both versions the parameter 7b = 270 MeV corresponds to the transition temperature in the pure 
YM theory. The third version, proposed by Fukushima, has only two parameters a and b and is 
inspired from a strong-coupling analysis 

^Fuku t>_ 

The parameter a determines the deconfinement transition in pure gauge theory, while b controls the 
mixing of the chiral and the deconfinement transition. Here, they are not fitted to the lattice data 
but also reproduce a first-order phase transition at 7b ~ 270 MeV in the pure YM theory. Details 
concerning the PQM model can be found in the forthcoming work [|I0|]. 

When dynamical quarks are present in the system the running coupling of QCD is modified 
by fermionic contributions. The size of this effect can be estimated within perturbation theory, see 

m 



e.g. As shown in Ref. [ ]12| ] this leads to an A^-modification of the expansion coefficients in 
the polynomial Polyakov-loop potential and can further be mapped onto an ./Vy -dependent Tq. In 
Tab. [I] one sees that the critical temperature 7b decreases for increasing Nf. 



N f 





1 


2 


2+1 


3 


T Q [MeV] 


270 


240 


208 


187 


178 



Table 1: The critical temperature 7b for Nf massless flavors according to [|12[|. The value for 2 
has been estimated by using HTL/HDL theory for a massive strange quark with m s ~ 150MeV. 



1 flavors 



2.2 Thermodynamic potential 

The thermodynamic potential of the PQM model is evaluated in mean-field approximation 
similar to [12, [l3|]. It splits into three contributions: the mesonic, U (a x ,a y ), the quark/antiquark, 
Q.q q , and the Polyakov loop contributions 

Q = U (a x ,a y )+a 9q (a x ,a y ^) ($,$) . (2.7) 

In general, the grand potential is a function of the temperature and three quark chemical potentials, 
one for each flavor. Here, we consider the isospin-symmetric case with two degenerated light quark 
masses. As a consequence, only two independent quark chemical potentials, the light pL q and the 
strange pL s , emerge. Furthermore, only two order parameters in the meson sector, the non-strange 
<j x and strange oy condensate are present. Note, the Polyakov loop variables are coupled to the 
fermionic part that also depends on the mesonic condensates via the quark masses. Finally, the 
temperature and quark chemical potential dependence of all four order parameters for the chiral 
and deconfinement transition are determined as solutions of the corresponding coupled equations 
of motion, i.e., 

dQ. _ da _ da _ da 

x y nun 

where min = [p x = (o x ) , G y = (o y ) ,<!> = ($) ,<!> = ( < l > )) labels the global minimum of the grand 
potential. 

In the following we have chosen a parameter set where both transitions coincide at /I = as 
suggested by lattice data which we use for comparison [Q]. This is the case for m a = 600 MeV and 
To = 270MeV where all three Polyakov-loop potentials yield coinciding transitions. 
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(a) Polyakov loop expectation value {<&) 
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(b) Subtracted condensate A c , s 



Figure 1: The Polyakov loop expectation value (<t>) (left panel) and the subtracted condensate A £/ . s (right 
panel) as a function of temperature for different Polyakov-loop potentials in comparison to lattice data for 
N T = 8 (symbols) [Q] . Solid lines correspond to larger pion and kaon masses as used in the lattice simulations, 
dashed lines to physical masses. 



3. QCD thermodynamics 



In Fig. |l(a)| the Polyakov loop expectation value (<&) is shown as a function of temperature 
for three different Polyakov loop potentials. For a proper comparison of the model results with 
the HotQCD [Q] and RBC-Bielefeld lattice data [|3j], we also adjust the pion and kaon masses 
accordingly to m n = 220 MeV and tuk = 503 MeV. Heavier meson masses yield also slightly 
heavier constituent quark masses. The increased meson masses shift the transitions also to higher 
temperatures. The strongest shift is seen in the non-strange chiral transition since the quark masses 
are largest in this sector and both transitions, the non-strange and deconfinement one, still coincide. 
The strange sector is almost unaffected. Solid lines in Fig. l(a)| correspond to heavier meson 
masses, used in lattice simulations and dashed lines to the physical masses. The values are slightly 
higher for the Fukushima potential than for the logarithmic potential which exhibits a sharper 
crossover. The influence of the higher meson masses is mild. The logarithmic potential comes 
closet to the lattice data in the broken phase, while the difference between the model results and 
the lattice is considerable in the restored phase for all potentials. 

For n = the Polyakov loop expectation value (<&) is available in model as well as in lattice 
calculations. This is in contrast to the chiral condensate (qq). A direct comparison of (qq) or 
(ss) with lattice simulations is difficult since unknown renormalization/normalization factors are 
involved. Therefore, a better quantity to compare with lattice simulations is the ratio 



(qq)(T) - (m q /m s )(ss)(T) 



(3.1) 



(qq)(0)-(m q /m s )(ss)(0) 

which involves the non-strange (qq) and strange (ss) condensates, correspondingly (a x ) and (o y ) 
in our model calculation. Since the bare quark masses rh q ,m s are not directly available, the ratios 



of the explicit symmetry breaking parameters of our model are used instead, cf. [15]. 



The results for A qjS are shown in Fig. 1(b). The model calculations are in better agreement than 
for the Polyakov loop expectation value. However, also here, the lattice data exhibit a smoother 
transition than the model results. The logarithmic potential generates the sharpest transition while 
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Figure 2: The normalized pressure as a function of temperature. The model calculations (PQM model with 
various Polyakov loop potentials and the QM model Jl4|]) are compared to lattice data similar as in Fig. 1. 



the Fukushima and polynomial potential are closer to the lattice results. The influence of the larger 
meson masses is stronger in the broken phase where the transition becomes smoother for larger 
meson masses. 

The pressure p(T,fi) = — £2 (T,n) is directly obtained from the thermodynamic potential with 
the normalization p(0,0) = 0. In Fig. || the (P)QM model pressures, normalized to the SB value of 



the PQM model, are compared to lattice data. As expected, the QM model [14] fails in describing 
the lattice data for all temperatures, while the PQM model results are in reasonable agreement with 
the data. The best agreement with the lattice data is achieved with Fukushima's potential but fails 
for temperatures above 1.57^. Physical meson (quark) masses (dashed lines) push the pressure to 
lower values, in particular around the transition. 

3.1 Taylor expansion and a novel algorithmic differentiation method 

As already mentioned one approach to overcome the sign problem in finite density QCD lattice 
simulations is based on an extrapolation from zero chemical potentials by a Taylor expansion. For 
the pressure the expansion in powers of (ju/r) at jx = reads 



'pA 



I>on 



n=0 



with c n (T) 



1 d"(p(T,n)/T A 



n\ d{}i IT)" 



(3.2) 



The coefficients c„(T) are accessible in QCD lattice simulations. Because of the CP-symmetry 
of the QCD partition function, Z(ju) = Z(— ju), only even coefficients contribute. These have been 
calculated by different lattice groups forNf = 2 [16, |T7| and Nf = 2 + 1 quark flavors [18], currently 
up to order n < 8. For the realistic 2+1 flavor scenario N r = 4, 6 lattices have been used and higher 



coefficients still have significant errors [18]. 

In principle, these coefficients can be obtained without any limitations in model calculations. 
Doing this analytically is a tedious and difficult process. Furthermore, standard numerical tech- 
niques like divided differences fail since the errors increase rapidly. In order to proceed we have 
developed a novel numerical technique based on algorithmic differentiation (AD). This technique 
allows the evaluation of higher derivatives to extremely high precision. In fact, it is essentially 
limited only by machine precision. Details of this method can be found in Ref. [^]. 
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0,98 1 1,02 0,98 1 1,02 0,98 1 1,02 

T/T X 

Figure 3: The Taylor coefficients C(, to cn in the PQM model with the logarithmic Polyakov-loop potential. 



In Fig. H we show the resulting coefficients C(, to C22 for the 2 + 1 flavor PQM model with 
the logarithmic Polyakov loop potential as a function of temperature. With increasing order n, the 
coefficients c n start to oscillate around the transition temperature T x . Inside a narrow temperature 
window around T x , i.e., in the interval 0.957^ >T> l.05T x the oscillation amplitude of the coef- 
ficients increases while they tend to zero far away of T x . This means that higher order coefficients 
become more and more important in the expansion even for fi/T < 1 and the convergence of the 
series becomes questionable. Nevertheless, the Taylor expansion might still be useful to gain some 
information on the existence and location of a possible critical endpoint (CEP) in the QCD phase 
diagram. For some recent estimates by lattice groups using only a few coefficients see [|TJ], [D|]. 
However, the validity of these results and the convergence of the series remain questionable. 

An often considered quantity to extract some information on the CEP is the quark number 
susceptibility % q = (d 2 Q.) /(djj, 2 ). Since % q diverges exactly at the CEP it is a suitable quantity to 
locate this point at least in model calculations. Its expansion can be expressed with the same Taylor 
coefficients c„ via 

^ = £„(„-iMr)(£)"-\ (3.3) 

n 

In Fig. || we compare % q obtained with the Taylor expansion of different orders with the evalu- 
ation of the model. For the chosen parameters the CEP is located at (T c ,^l c ) ~ (185, 167)MeV, i.e., 
at (jx c /T c ) ~ 0.9 (middle panel). The other two panels show % q for the ratios (ju/r) = 0.8 and 1.0, 
respectively. As expected, % q diverges at T c ~ 0.9T X and is finite for larger or smaller ratios (jx/T) 
near the critical one (lines labeled with 'PQM'). Even the 12 th order expansion neither reproduces 
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(c) n = 4,8, 12, (£) = 1.0 



Figure 4: The quark-number susceptibility % q /T 2 for different ratios of \ijT and different orders of the 
Taylor expansion. The black line, labeled with 'PQM', is the model calculation. 



this structure nor shows any significant remnant of the divergence at the CEP. In the chirally broken 
phase, i.e., on the left side of the peak, the agreement becomes better with increasing orders of the 
Taylor expansion. The peak structure which is visible in the Taylor expansion around T ~ T x is a 
pure truncation artifact and a signal for its breakdown near the transition. On the other hand, in the 
chirally restored phase at high temperatures, the expansion again reproduces the model result, even 
for jtt/r = 1. This is in agreement with the high-temperature limit of the coefficients c n since only 
the coefficients with n < 4 have a finite Stefan-Boltzmann value. 

For a deeper understanding of the breakdown of the Taylor expansion it is instructive to con- 
sider its convergence radius. It can be obtained from the definition 



r = lim rm = lim 



Cln 



C2n+2 



1/2 



(3.4) 



It is not known how well the radius r is estimated by r„ for a finite value n. 

In Fig. H we show the convergence radius for different orders n in the phase diagram together 
with the chiral and deconfinement phase boundaries. For temperatures close to T x the oscillations 
in the coefficients cause oscillations in r n . For smaller temperature, i.e., approaching the CEP from 
above, the CEP lies inside the convergence region for lower n values. For larger n values r n (T c ) 
approaches the CEP and points to jj, c once T c is known. In the first order phase transition region, 
i.e., for T < T c the interpretation of the convergence radius r n becomes misleading. At a first order 
transition a new global minimum appears in the grand potential and the order parameter jumps 
to the new minimum. But the Taylor expansion is still performed around the /i = minimum 
and cannot capture the jump. This means that its convergence radius is only valid for this pL = 
minimum. Therefore the results obtained with the expansion are only valid up to the first order 
phase boundary. 

Once T c is known, the convergence radius might provide an estimate for the critical pL c and 
to locate the CEP. In Jl7| ] the critical temperature was estimated by using the sign of the Taylor 
coefficients. It contains information on the location of the CEP in the complex chemical potential 
plane [|l^]. For temperatures below the critical one all coefficients should be positive and the first 
root in the coefficients for T < T x should determine T c in the limit n — > oo. Since the location of 
the CEP is precisely known in the model calculation we can verify this idea. From the coefficients, 
shown in Fig. we observe that the first roots of the coefficients, starting at T ~ T x for n = 6, 
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Figure 5: The convergence radii r„ for different orders n of the expansion in the PQM model with the 
logarithmic Polyakov loop potential (7b = 270 Me V). Also shown are the phase boundaries for the chiral 
transition (red line; dashed: crossover, solid: first order) and the deconfinement crossover transition (yellow 
line). The black dot indicates the CEP. 



stay far away from the critical temperature T c ~ 0.9T X . From this one might conclude, that the 
Taylor expansion yields a too large value for the critical temperature. A more detailed study in this 
direction is underway [p0|]. 

4. Summary 

In this talk we have presented recent work, which extends a two flavor Polyakov-quark-meson 
model to Nf = (2 + 1) quark flavors. The bulk thermodynamic of this model in mean-field approx- 
imation is compared to recent (2 + l)-flavour QCD lattice simulations. The larger quark masses on 
the lattice are also considered in the model calculations. For the Polyakov-loop potential we exam- 
ine three different ansatze which all reproduce a first-order phase transition in the pure gauge sector 
at the critical temperature T c = 270 MeV. For the used parameter sets of the Polyakov-loop poten- 
tials the chiral and the deconfinement transition coincide at vanishing chemical potential. With the 
Polyakov-loop a very good agreement of the QCD equation of state in particular in the transition 
region up to temperature of T ~ 1-57^ is achieved in contrast to a pure three flavor quark-meson 
model. 

Furthermore, we analyze the finite chemical potential expansion of the quark-number suscep- 
tibility in a Taylor series around vanishing chemical potential. By means of a novel algorithmic 
differentiation technique, we calculate the Taylor coefficients up to 24 th order in the PQM model 
for the first time. The knowledge of these higher Taylor coefficients allows the systematic study of 
convergence properties of the series. The peak in the quark-number susceptibility % q , seen in the 
Taylor expansion, is a pure artifact and related to the breakdown of the expansion. In particular, 
the divergence of % q at the critical endpoint in the phase diagram cannot be captured by the Taylor 
series. Nevertheless, away from the phase transition the Taylor expansion converges rapidly and 
coincides with the model evaluation even for \ijT = 1. In the first-order region discontinuities in 
the order parameter emerge which cannot be captured by the Taylor series either. However, our 
analysis of the convergence radius yields that the expansion might work even for jj./T > 1 if one 
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stays in the broken phase. But in lattice simulations this transition boundary is not known. Thus, in 
order to make statements on the possible location of a critical endpoint in the phase diagram with 
the Taylor expansion technique our studies suggest that higher order coefficients with n > 8, more 
as currently available on the lattice, are definitely required. 
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